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PARA 


A  COMPUTER  SIMULATION  CODE  FOR  PLASMA 


DRIVEN  ELECTROMAGNETIC  LAUNCHERS 


1 .  INTRODUCTION 


This  report  presents  the  mathematical  details  of  PARA,  a  computer 
code  for  simulation  of  plasma  ^rmature  r_ail  accelerators.  Outlines  of  the 
code  have  been  given  previously  by  Thio  (1980  [11,  1981  [21  ,  1982  [3]). 


Railgun  performance  simulation  codes  have  proliferated  in  the  past 
few  years  [4-71  owing  to  the  current  military  and  industrial  interests  in 
railguns  and  the  role  they  play  in  railgun  design  and  system  studies. 
However,  to  the  best  knowledge  of  the  author,  none  of  the  simulation  codes, 
other  than  PARA,  has  the  capability  of  modeling  and  predicting  the  detailed 
physical  properties  of  the  plasma  armature  as  well  as  performance.  The 
earlier  codes  [5,6]  assume  a  constant  value  of  200  v  as  the  voltage  drop 
across  the  plasma  armature.  More  recent  codes  [7]  allow  for  a  time-varying 
voltage  drop  by  assuming  a  fixed  resistance  for  the  plasma.  Other  details 
such  as  temperature,  pressure  and  length  of  the  plasma  are  not  calculated. 


The  first  serious  attempt  at  studying  the  behaviour  of  the  plasma 
armature  was  made  by  McNab  [8].  Following  that,  two  independent  studies  were 
launched  almost  concurrently:  One  by  Powell  and  Batteh  [9]  at  the  Ballistics 
Research  Laboratory,  Aberdeen  (1979);  the  other  by  Thio  [1-3]  at  the 
Materials  Research  Laboratories,  Maribyrnong  (Australia)  where  a  parallel 
integrated  experimental  study  of  the  plasma  armature  was  also  undertaken 
(1980). 


ttie  work  of  McNab  was  essentially  a  steady  state  analysis.  Taking 
McNab's  work  as  a  first  approximation,  Powell  and  Batteh  [9]  extended  the 
analysis  but  still  retained  the  steady-state  feature.  Specifically,  Powell 
and  Batteh  assumed  constant  current  and  a  plasma  column  in  thermodynamic 
equilibrium.  Later,  Powell  [10]  extended  the  analysis  to  2-D.  No  attempt 
was  made  to  incorporate  the  model  in  a  simulation  code.  Thus  no  effects 
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associated  with  the  non-linearity  of  the  arc  as  a  circuit  component  of  the 
driving  railgun  circuit  and  no  temporal  variation  of  the  plasma  properties 
were  studied. 


PARA  represents  the  first  attempt  at  modeling  some  of  the  time- 
varying  properties  of  the  plasma  by  using  a  quasi -static  approach.  More, 
assumptions  were  made  for  the  plasma  model  in  PARA  than  by  Powell  in  the 
interest  of  mathematical  simplification.  However,  in  PARA,  the  plasma  model 
is  embedded  in  a  comprehensive  computer  code  for  the  simulation  of  railgun 
operation.  This  enables  one  to  follow  the  dynamical  history  of  the  plasma  in 
a  typical  rail-launcher  firing  and  study  the  time  development  of  the  various 
plasma  parameters  including •  voltage  drop,  volume,  ten?>erature  and  pressure  of 
the  plasma.  Knowledge  of  these  plasma  properties  is  particularly  important 
in  a  systematic  treatment  of  the  material  problems  associated  with  the 
eventual  realization  of  plasma-driven  electromagnetic  launcher  technology. 


•  The  circuit  equations  used  in  PARA  are  introduced  in  Section  2  where 
the  major  circuit  components  are  identified.  In  Section  3,  the  details  of 
the  plasma  model  are  presented.  For  convenience  of  reference,  all  pertinent 
equations  used  in  PARA  are  assembled  in  Section  4,  and  the  broad  structure  of 
the  computer  program  is  outlined  in  Section  5.  The  concluding  section 
contains  a  critique  of  the  present  version  of  PARA  and  some  recommendations 
for  future  development  of  the  code. 


2.  CIRCUIT  EQUATIONS 


2.1  The  Circuit 


The  equivalent  circuit  used  for  the  railgun  is  shown  in  Figure  1. 
The  power  supply  consists  of  a  capacitor  bank  C,  which  is  connected  in  series 
with  an  inductor  I*0  via  a  switch  S1  which  may  be  an  ignitron  or  a  spark 
gap.  The  switch  S,  is  the  main  switch  which  turns  on  the  main  discharge. 

The  capacitor  bank  is  also  equipped  with  a  crow-bar  switch  Sj,  which  in  the 
code  can  be  sat  to  turn  on  at  any  particular  voltage  across  the  capacitor. 

In  practice  depending  on  the  actual  switch  used,  there  are  constraints  on  the 
bias  voltage  which  Is  needed  to  activate  it. 

The  rails  are  connected  in  series  with  the  inductor  1>Q.  The 
reactance  in  the  rails  is  lumped  as  a  variable  inductance  Lr  and  a  variable 
resistance  Rr.  The  plasma  armature  manifests  itself  as  a  variable  resistor 
Rp  between  the  rails  with  Vg  g  and  VE  g  the  electrode  potential  drop  at  each 
rail.  Similarly,  for  the  plasma  switches  S,  and  S2,  a  resistance  (Rg  ,  or 
Rs  j'  *nd  electrode  potential  drops  (VE  VE  2,  3  *nd  vs  4)  ar*  allowed 

for.  In  the  case  of  the  switches,  the' resistance  is  assumed  to  be  constant. 
Any  stray  resistance  associated  with  the  bus-bars  and  with  any  electrical 
connections  in  the  capacitor  crow-bar  switch  circuit  are  lumped  with  the 
resistances  Rg  1  and  Rg  2.  Stray  resistance  across  contact  joints  in  the 
main  railgun  circuit  are  lumped  as  a  fixed  resistor  Rg  3.  The  resistance 
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of  the  main  inductor  L0,  hcwever,  is  allowed  to  vary  with  time,  in  an  attempt 
to  take  account  of  the  skin  effect.  Any  stray  inductance  in  the  main  driving 
circuit  is  also  lumped  with  t,  . 


In  the  computation,  the  coupling  between  the  capacitor-crow-bar 
circuit  and  the  main  driving  circuit  after  the  crow-bar  Sj  is  activated  is 
completely  ignored.  Thus,  after  S-  is  turned  on,  computation  proceeds  as  if 
the  capacitor  bank  has  been  physically  removed  from  the  circuit,  and  the 
circuit  consists  of  inductances  and  resistances  only. 

The  electrode  drops,  V_  , ,  ,,  V_  -  and  V_  ,  are  lumped  as  a 

,  _  .  1  t  *,  1L  f  D  Cif  O 

single  voltage  Vg  0,  given  by 


E,  B 


V  V  +  V  +  V 

E,1  E,2  E,5  E,6 


(2.1  ) 


Similarly,  after  crow-bar. 


E ,  A 


V  +  V  +  V  +  V 
E,3  E,4  E,S  E,6 


(2.2) 


The  stray  resistances  are  lumped  thus, 


R„  _  »  R_  ,  +  R„  , 

S, B  S,  1  S ,  3 


(2.3) 


before  crow-bar,  and 


RS,A  “  RS,2  +  RS,3 


(2.4) 


after  crow-bar.  In  the  current  version  of  PARA,  it  is  assumed  that 


E,  A 


v_  „  i-e.  v.  ,  +  v  ,  >  v  ,  +  v„  „ 

E ,  B  E ,  1  E,  2  E ,  3  E ,  4 


and 


(2.5) 


R„  „  “  R_  _  i.e.  R„  ,  ■  R„  . 

S , A  S,B  S,1  S ,  2 


2.2  The  Equations 


With  switch  S1  closed  and  Sj  open,  the  basic  equations  governing  the 
temporal  behaviour  of  the  circuit  are: 
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(2.6) 


dl 

dt 


VC 


I(RS  +  Rb  *  Rr 

L  +  L'Z 
o 


+  R  +  L  '  V 

££21 


dV 

_ c 

dt 


(2.7) 


where  I  is  the  current  in  the  rails  and  through  the  plasma  armature  and  Vc  is 
the  voltage  across  the  capacitor;  VE  is  VE(A  or  VE  B  and  Rg  is  Rg  A  or  Rg  B; 
Rb»  Rr  and  Rp  are  variable  resistances  of  the  inductor  LQ,  the  rails,  and  the 
plasma  column  respectively.  L'  is  the  inductance  per  unit  length  of  the 
rails,  z  is  the  position  of  the  back-errd  of  the  projectile  and  V  ,  i3  its 
velocity.  P  3 


After  the  switch  S2  closes,  the  circuit  equations  (2.6)  and  (2.7) 
reduce  to  a  single  equation 


dl 

dt 


E,  B 


I(R 


5 ,  B 


+  R,_  +  R  +  R  +  L 


V  .  ) 

Pr°J- 


+  L'z 


(2.8) 


For  calculating  the  resistances  Rb  and  Rf,  the  conductors  for  the  inductor, 
the  bus-bars  and  the  rails  are  assumed  to  have  rectangular  cross-sections  or 
cross-sections  which  can  be  derived  from  a  small  amount  of  conformal 
deformation.  A  skin  depth  5  is  assumed  for  these  conductors 
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(2.9) 


which  is  associated  with  a  step  current  pulse  [11,121,  where  n  is  the 
resistivity,  u  permeability  and  t  is  time.  The  usual  simple  expression  for 
resistance 


c 


is  used  for  calculating  Rv  and  Rr,  where  l  is  the  length,  wc  the  width  and  d 
is  either  the  thickness  of  the  conductor  or  twice  the  skin  depth  26,  whichever 
is  smaller.  Velocity-skin  effect  for  the  rails  is  not  considered  in  the 
present  version  of  PARA. 


The  yet  undefined  quantities  in  the  governing  equations  (2.6)  to 
(2.8)  are  the  position  z  and  the  velocity  v  ,  of  the  projectile,  and  the 
resistance  Rp  of  the  plasma  column,  which  are  treated  in  the  next  section. 
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3.  PLASMA  MODEL 


3. 1  General  Assumptions 


The  following  assumptions  are  made: 

(1)  The  conditions  of  LTE  (Local  Thermodynamic  Etjui librium ) ,  which 
mean  that,  (a)  the  various  particle  species,  electrons,  ions  and  neutral 
particles  all  share  the  same  temperature  at  any  particular  instant  of  time: 
and  (b)  Saha's  equation  may  be  used  for  calculating  the  degree  of  ionization. 

(2)  The  Debye  length  is  sufficiently  small  for  a  one-fluid  MHD 
approach  to  be  appropriate. 

(3)  The  electrical  conductivity  of  the  plasma  is  calculated  with 
Spitzer-type  expressions.  These  expressions  assume  a  high  degree  of 
ionization  and  ignore  the  presence  of  a  magnetic  field. 

(4)  Effects  of  viscosity  are  neglected. 

(5)  For  field  quantitfi.es,  only  variations  in  the  direction  parallel 
to  the  rails  are  considered. 

(6)  Current  density  and  temperature  are  assumed  uniform  throughout 
the  plasma. 

(7)  The  mass  of  the  plasma  is  assumed  to  remain  constant. 

(8)  The  inertia  force  in  the  plasma  is  assumed  to  be  very  much 
smaller  than  the  Lorentz  force  and  the  pressure-gradient  force. 

(9)  The  plasma  is  assumed  to  be  singly  ionized  from  a  single  atomic 

species. 


3. 2  Geometry 

The  geometry  is  as  shown  in  Figure  2.  It  is  assumed  that  the 
plasma  is  confined  to  a  rectangular  region  between  z  *  za  and  z  »  z^,  the 
current  flows  through  the  plasma  in  the  negative  y-direction  with  the  current 
density  vector  ■  -  j  e^,,  the  magnetic  induction  B_  points  in  the  positive 
x-direction  (into  the  plane  of  paper),  and  the  plasma -projectile  travels  in 
the  positive  z -direction  with  the  velocity  V_  »  vproj  ®x*  field 

quantities  are  assumed  to  vary  only  in  the  z -direction. 

3. 3  Magnetoplasmadynamic  Equations 

Following  Ferraro  and  Plumpton  [13J,  the  magnetohydrodynamic  (MHD) 
equation  for  momentum  transport  is 

P^  +  7p  ■  jxB+4  pv7  (V  .  v)  +  pvV2  v  (3. 1  ) 

at  “  j  —  — 

where  p  is  mass  density  of  the  plasma,  _v  the  velocity,  p  is  the  pressure  which 
includes  both  the  gas  and  radiation  pressure,  v  is  the  kinematic  coefficient 
of  viscosity. 
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Viscous  effects  are  ignored  (v  »  0)  and  it  is  assumed  that 


P  ~  <<  Vp  (3.1a) 


and  so,  P  —  «  i  x  8  (3.1b) 

This  is  approximately  the  same  as  saying  that  the  inertial  force  acting  on  the 
plasma  is  small  compared  to  the  overall  Lorentz  force  which  accelerates  both 
the  plasma  and  the  projectile.  With  these  assumptions,  equation  (3.1) 
reduces  to 


Vp  »  ”  «  2.*  ! 


(3.2) 


Use  is  made  of  Ampere's  law  from  Maxwell's  equations  to  relate  the 
current  density  to  the  magnetic  field 


V  x  H  -  J 


(3.3) 


In  the  coordinate  system  of  Figure  2,  equations  (3.2)  and  (3.3)  have  the  form 


H  -  u  1“ 

(3.4) 

3H 

(3.5) 

37  "  '5 

It  is  then  assumed  that  the  current  density  is  uniform  in  spa'-e 
throughout  the  plasma  so  that  equations  (3,4)  and  (3,5)  can  be  integrated 
directly  to  give  the  magnetic  field  within  the  plasma 


H(z,t)  -  H(za,t)  -  j(t)  (z  -  z#)  (3.6) 

and  the  pressure  distribution  in  the  plasma 

p(z,t)  *  p(z  ,t)  +  u  j(t)  {H(z  ,t)  (z  -  z  )  -  (z  -  z  )2}  (3.7) 

The  pressure  in  the  plasma  is  a  superposition  of  two  terms:  a  term  p(za,t), 
corresponding  to  the  pressure  at  the  rear  boundary  za,  and  a  term  (2nd  term  in 
(3.7))  associated  with  the  current  distributed  in  the  plasma. 


p 


RT(1  +  a) 


(3.9) 


where  a  the  degree  of  ionization,  is  used  with  expression  (3.7),  to  evaluate 
the  mass  integral  to  give 


*  “t1  '"<Ve>  -  ¥■'! 


(3.  10) 


where  l  is  the  length  of  the  plasma  armature.  In  the  above  integration, 
temperature  T  has  been  assumed  to  be  uniform  in  space.  Despite  the 
uniformity  of  temperature,  the  degree  of  ionization  is  in  general  not  uniform 
according  to  Saha's  equation  if  the  density  and  pressure  also  vary.  However, 

in  obtaining  the  result  (3.10),  the  degree  of  ionization  a  is  assumed  to  be 
uniform  in  space,  again  in  the  interest  of  mathematical  simplicity. 
Rearranging  expression  (3.10),  the  following  expression  for  the  pressure  at 
the  rear  boundary  is  obtained; 


P(za.t) 


m  RT(t) (1  +  a(t) ) 
P _ 


If  the  plasma  is  pushing  against  a  closed  breech,  the  pressure  at 
the  rear  boundary  p(za,t)  is  non-zero,  giving  rise  to  an  additional  force  on  - 
the  projectile  due  to  the  thermal  expansion  of  the  plasma.  the  term  p(za,t) 
is  calculated  as  part  of  the  computation  and  in  this  sense  the  'explosive' 
propulsion  force  due  to  the  expansion  of  the  plasma  is  accounted  for  in  PARA. 
The  force  of  explosion  associated  with  the  transition  between  solid  and  vapour 
during  the  fuzing  of  a  metallic  foil  to  generate  the  plasma  is  not  calculated 
in  PARA.  The  energy  associated  with  this  phase  of  explosion  is  however  very 
small  and  the  purturbation  caused  by  it  to  the  overall  dynamics  of  the  plasma- 
projectile  should  be  negligible  even  in  the  low  energy  program  RAPID  at  MRL 
(2,3). 

When  the  plasma  armature  is  completely  confined  by  magnetic  pressure 
and  not  in  contact  with  the  breech,  the  pressure  p(za,t)  at  the  rear  boundary 
is  zero.  In  this  case,  expression  (3.11)  yields  a  relationship  between  the 
length  of  the  plasma,  the  temperature,  degree  of  ionization,  the  mass  of  the 
plasma,  and  the  current  in  the  plasma. 
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if/v.V-v.V/.  The  latter  is  approximately  the  case  if  the  variation 
in  the  velocity  field  is  small,  and  this  is  assumed  to  be  the  case.  It  is 

also  assumed  that  the  plasma  expands  and  contracts  uniformly,  that  is  div  v_  is 
spatially  uniform,  so  that, 


div  v  =  V/V 


(3. 16) 


where  V  is  the  volume  of  the  plasma.  With  these  approximations,  (3.14)  is 
reduced  to 


d£  _  dQ  _  dV 

dt  dt  ^av  dt 

where  U  is  the  total  internal  energy  of  the  plasma, 


(3.17) 


CJ  =  /v  u  dV  . 

the  quantity  pav  is  the  volume  average  of  the  pressure  field. 


P 

av 


_1_ 

v 


/v  p  dV 


=  p ( z  ,t)  +  [H (z  ,t)  -  ii]  (3.18) 

A  *  3  w 


upon  substituting  expression  (3.7)  for  the  pressure  field  within  the  plasma; 
or  using  the  equation  of  state  (3.9)  and  using  the  same  approximation  as  in 
evaluating  the  mass  integral,  Pav  can  be  derived  to  be. 


P 


av 


(1  +  a)  m  RT 

_ E _ 

V 


(3.19) 


The  quantity  Q  is  defined  by 


dQ 

dt 


d_ 

dt 


/v  e  dV 


/  V 
V 


F  dV 


i2r 


3V  —  *  - 


ds 


(3.20) 


where  the  surface  integral  gives  the  energy  loss  by  radiative  flux  through  the 
boundary  of  the  plasma.  This  is  computed  assuming  the  plasma  radiates  like  a 
b  lack -body : 

F  .  n  dS  =  2 (wh  +  hi  +  Iw)  a  T4 
3  V  s 
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where  g  is  Stefan's  constant  for  black-body  radiation,  T  is  the  plasma 
temperature,  h  is  the  height  of  the  bore  and  w  is  the  width  of  the  bore 
(separation  between  the  rails). 

Equation  (3.17)  can  now  be  treated  unambiguously  in  terms  of  the 
total  internal  energy  of  the  plasma.  For  this  the  following  form  for  the 
internal  energy  is  assumed; 


0 


3  4 

—  m  ST  1  +  a  +  a  m  e  +  a  VT 
2  P  pi 


(3.21) 


which  is  appropriate  for  a  singly  ionised,  monoatomic  species  with  e^  as  the 
ionization  energy  per  unit  mass,  and  a  is  the  Stefan-Boltzman  constant.  The 
corresponding  Saha  equation  for  ionization  is 


exp 


,  q  V  . 
{*  } 


(3.22) 


where  Mr  is  the  atomic  weight,  NQ  Avogadro's  number,  m0  electron  mass,  k 
Boltzman  constant,  h  Planck's  constant,  qg  electron  charge,  and  the 
ionization  potential. 

the  total  time  rate  of  change  of  the  internal  energy  may  be 
expressed  as 


dU  dT 

dt  “  (3Tlv  dt  + 


d V 
dt 


Hie  partial  differential  coefficients  may  be  calculated  from  (3.21)  and 
(3.22),  and  using  E^/R  *  VJA,  to  give 


.  /  _3  g(  1  -  op  /  3  qe  VI  \  -  ?r  \ 

a)  \2  *  (2  -  a)d  +  a)  (2  +  kT  /  +  12  p^  I 

iiHl  -  ad  -  a)  /3  qe  VI  \ 

\  3V  L  Pr  ( 2  -  a)  ( l  +  a)  (  2  kT  f  pav 


mp  R(  1  + 


(3.23) 


(3.24) 
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where  p  *  —  aT  is  the  radiation  pressure.  It  is  worth  noting  that,  by 
definition,  the  specific  heat  per  unit  mass  for  the  plasma  at  constant  volume 
is 
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which  i3  expression  (3.23)  apart  from  the  factor  m  .  Expression  (3.23)  is 
thus  an  expression  for  the  specific  heat  cv  for  a  plasma  which  takes  into 
account  the  energy  for  single  ionization  and  radiation  pressure  and  whose 
equation  of  state  is  that  for  an  ideal  gas  modified  for  ionization.  In  term 
of  c  ,  the  heat  development  equation  is. 


dT 

mp  cv  dt 


d£  dV 

dt  E1  pav  dt 


(3.25) 


where 


a(i  -  oO 

+  ( 2  -  a)  ( 1  +  a) 


%  V 

kT 


+ 


3  P. 


(3.26) 


Comparing  with  (3.12),  it  can  be  seen  that  an  unambiguous  interpretation  of 
the  specific  heat  c,  and  the  pressure  p  with  the  factor  f1  as  a  'correction' 
factor  has  been  derived. 


3.5  Current  Density  and  Boundary  Magnetic  Field 


Before  expression  (3.7)  can  be  applied  to  obtain  the  gas  pressure  at 
the  base  of  the  projectile,  expressions  for  the  current  density  j(t)  and  the 
magnetic  field  at  the  rear  boundary  H(za,t)  are  required.  By  assumption  of  a 
uniform  current  density  within  the  plasma 


j(t) 


(3.27) 


where  I  is  the  total  current,  4  the  length  of  the  plasma  armature,  h  the 
height  of  bore,  and  (1  -  y  )  the  fraction  of  current  which  is  lost  through 
arcing  ahead  of  the  projectile;  the  factor  y 2  is  normally  set  to  1  in  the 
running  of  the  code. 

For  calculating  the  magnetic  field  H(za,t)  it  the  rear  boundary,  a 
simple  but  rigorous  approach  is  to  assume  the  case  of  infinitely  high  rails. 

In  that  case,  it  can  be  shown  that  the  magnetic  field  in  the  region  between 
the  rails  sufficiently  far  downstream  from  the  plasma  is 

H(z  ,t)  -  £  (3.28) 

3  n 

This  is  the  approach  used  by  Powell  and  Batteh  [9],  It  does,  however,  tend 
to  overestimate  the  average  field  between  the  rails. 

A  less  rigorous,  but  perhaps  more  realistic,  approach  is  to 
calculate  the  field  in  the  region  between  the  rails  due  to  thin  sheets  of 


current  flowing  in  the  inner  surface  of  the  rails  of  finite  height  and  average 
the  field  over  the  bore.  This  approach  is  adopted  in  the  current  version  of 
PAHA,  but  instead  of  averaging  the  field  over  the  bore,  it  is  only  averaged 
along  the  mid-line  between  the  rails. 

The  algebraic  result  is 


H(ta,t)  -  {r1(h/2w)  +  Y2/2}  I/h  (3.29) 

where  r,  is  a  function  given  by 

y ,  (8 )  -  7  itan_1  8  +  j  Sin  (i  +  8"2)} 


3.6  Propulsive  Force 

From  (3.7),  the  pressure  at  the  base  of  the  projectile  is  given  by 
p(zb,t)  -  p(za,t)  +  Uj (t)  { H(z  , t )£  (3.30) 

Substituting  (3.27)  and  (3.29)  into  (3,30) 

,2 


p(zb,t)  ■  p  (z  a ,  t )  +  u72Y1(h/2w)  (I/h) 


(3.31) 


The  accelerating  force  acting  on  the  projectile  is  then  given  by 


p(z.  ,t)  wh 
D 


p(za>t)  wh  +  (liW/h)  y2y1  (h/2w)  X 


(3.32) 


This  expression  is  used  in  the  elementary  equation 


dfz 

dt2 


F 

_E 

m 


to  calculate  the  position  z  of  the  projectile. 


(3.33) 


3. 7  Plasma  Resistance 


The  resistance  of  the  plasma  armature  is  calculated  as: 
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n  w 


(3.  34) 


R 

P 


hi 


where  n  is  the  resistivity  of  the  plasma.  Spitzer's  expression 
strongl?  ionized  gas  in  the  absence  of  a  magnetic  field  and  ionic 
unity  (2  »  1)  is  used  to  calculate  the  resistivity: 


[Ml  for  a 
charge  of 


n 

P 


65.3  in  A 

t3/2 


13.  m 


(3.35) 


where 


1.238  x  107  T3/2/n  1/2 
e 


and  is  the  electron  density.  Though  the  expression  for  resistivity  is 
dependent  on  the  electron  density  which  varies  as  the  mass  density  p,  the 
dependence  is  very  weak  and  an  average  value  for  the  electron  number  density 
is  used; 


Total  number  of  electrons  in  the  plasma 
e  Volume  of  plasma 

a  N  m 

-  - 2 — E. 

M  V 
r 


(3.36) 


where  Nq  is  Avogadro's  number,  the  atomic  weight  of  the  atomic  species  and 
V  the  volume  of  the  plasma. 


4.  ASSOIBLY  OF  EQUATIONS 


Por  convenient  reference,  this  section  is  a  collection  of  the  basic 
set  of  equations  used  by  PARA.  Equation  numbers  from  the  text  are  shown 
wherever  applicable. 


V  -  I(R  +  R.  +  R  +  R  +  L’V 
E _ s  b  r  p _ pro]  ) 

L  ♦  L'z 
o 


(2.7) 


dV 

c 

dt 


1 

C 


(4.  1  ) 


(4.2) 


(3.32) 


dV 


££2l 

dt 


(p(za,t)  wh 


♦  (liw/h) 


Vi 


(h/2w) 


I2}  /(m_ 


+  M) 


(4.3) 


(4.4) 


1  3 


dz 

dt 


V 

pro] 
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I 


(3.34) ,  R 

(3.35)  F 


w  65.3  „  ,,  ,  .7  _3/2  ,  1/2. 

kT  ~3/2  ln  °*238  x  10  T  /"e  > 


(3.36)  n 


aN  m 
.  —  2— £ 
e  MV 
r 


(3.29)  L’  «  (2yw/h>  Y1(h/2w> 


(3.29)  y,(S)  -  i{tan'16  +  ~  6  ln  (1  +6-/)} 


(3.23)  c  -  R(1  +  a )  (|  +  _ - _  (A  + 

v  1 2  ( 2  -  a )  ( 1  +  a )  2  kT 


q  VT  2  .  4 

X,  .  4  aT 


~)  + 


'  av 


(3.26)  f,  -  1  + 


a(1  -  a)  r£ 

(2  -  a) (1  +  a)  ‘ 2 


3  ♦ 


kT 


aT 


m  RT(t)d  +  a(t)  ) 
(3.11)  p(za,t)  «  max  { 0,  -£ - - 


-  -U3-( *1  1  ~  x  [H(z  ,t)  -  -!—■ C— ’  ]} 


(3.29)  H(za,t)  -  {y ,  (h/2w)  +Y2/2}  I/h 


(3.22)  a2(1  -  a)"1  -  (M  /N  )  ( 2**  kT/h2)  3/2  (V/m  )  exp  t-q  VtAT] 

r  o  e  pel 


The  physical  meaninqs  of  the  symbols  are  given  below: 

X  Current  in  the  railgun  circuit 

Vc  Voltage  across  the  capacitor  bank 

VE  Total  electrode  potential  drop  in  the  circuit 

Rg  Stray  resistance,  constant  in  time 

Rjj  Time -dependent  resistance  of  the  inductor 

Rf  Time -dependent  resistance  of  the  rails 


(4.14) 


(4.15) 


(4.16) 


(4.17) 


(4.18) 


(4.19) 


(4.20) 


(4.21) 


(4.22) 


Rp 

L' 

Lo 

z 

c 

V 

pro] 

w 

h 

u 


Qr 


®rad 

Pav 


Time-dependant  resistance  of  the  plasma 
Inductance  per  unit  length  of  the  rails 
Inductance  of  the  storage  or  pulse-shaping  inductor 
Position  of  projectile 
Capacitance  of  the  capacitor  bank 
Projectile  velocity 

Pressure  at  the  rear  boundary  of  the  plasma 
Rail  separation 

Height  of  that  portion  of  rail  exposed  to  the  plasma 
Permeability  of  free  space 
A  geometric  function 

Fraction  of  current  lost  through  arcing  ahead  of  the  projectile. 
Normally  assume  »  1. 

Tesperature  of  plasma 

Stefan  constant  for  black -body  radiation 
Heat  dissipated  in  plasma  through  ohmic  heating 
Heat  dissipated  in  rails  through  ohmic  heating 
Energy  dissipated  by  electrode  potential  drop 
Energy  radiated  by  plasma 
Mean  pressure  within  the  plasma 
Specific  heat  of  the  plasma 
A  thermodynamic  function 

Resistivity  of  the  material  of  which  the  inductor  is  made 

Rail  resistivity 

Conductor  length  of  inductor 

Effective  thickness  of  inductor 

Effective  width  of  inductor 
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4^  Length  of  rail  up  to  the  position  of  projectile  from  the  breech 

wf  Total  rail  height 

ne  Electron  number  density 

a  Degree  of  ionization 

qe  Electronic  charge 

Vj.  Ionization  potential  for  the  atomic  species  making  up  the  plasma 

k  Boltzman  gas  constant 

a  Stefan-Boltzman  constant 

nip  Total  mass  of  plasma 

M  Mass  of  projectile 

Mr  Molecular  or  atomic  weight  of  the  ionic  species  in  the  plasma 

Nq  Avogadro  's  number 

me  Mass  of  electron 

h  Planck  's  constant 

S  Skin  depth 

Conusant:  the  dimensions  t  ,  d,  ,  wfa  are  specifically  designed 

for  the  case  where  the  main  inductor  is  a  coil  made 
up  of  flat  metallic  strips. 


5.  STRUCTURE  OF  THE  PROGRAM 


The  overlay  structure  of  the  program  is  shown  in  Figure  3.  The 
names  of  the  subprogram  units  are  displayed  within  rectangular  enclosures. 

The  paths  of  subprogram  calls  are  indicated  by  solid  lines  and  the  arrows 
indicate  the  directions  of  call.  A  cross  +  is  a  joint  in  the  path  whereas  a 
semicircle  ft  indicates  a  bypass.  Thus,  the  main  program  is  named  PARA  and  it 
calls  upon  four  subroutines:  SETUP,  NORMAL,  EXTRAC  and  SGUN.  Subroutine 
SETUP  handles  the  input  of  data  and  sets  up  the  data  array  CONDAT.  NORMAL 
normalizes  the  data  in  CONDAT  to  dimensionless  quantities.  EXTRAC  extracts 
the  data  required  for  the  first  gun  segment  (PARA  can  handle  guns  of  more  than 
one  segment)  and  stores  them  in  the  array  GUNDAT,  ready  for  use  by  SGUN. 

SGUN  calls  upon  INTER?  to  interpret  the  data  in  GUNDAT,  CAPSEG  to  simulate  a 
capacitively  driven  gun  segment,  and  INDSEG  to  simulate  an  inductively  driven 
gun  segment.  CAPSEG  calls  upon  CAPAC  for  a  definition  of  the  system  of 
equations  which  are  essentially  equations  (4.1 )— (4. 9)  and  calls  upon  RESULT  to 
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channel  the  results  of  the  computation.  Similarly  INDSEG  calls  upon  INDUCT  to 
obtain  the  equations  appropriate  for  an  inductively  driven  segment  which  are 
essentially  the  same  equations  (4.1) -(4. 9)  with  V.  =  0  and  dVc/dt  =  0.  It 
also  calls  upon  RESULTS  to  output  the  results  of  the  simulation.  CAPSEG  also 
calls  upon  INDSEG  after  the  capacitor  bank  is  crow-barred.  Both  INDUCT  and 
CAPAC  make  use  of  VETROD,  VRAIL  and  PLASMA  to  calculate  the  relevant 
quantities  from  the  list  of  definitions  (4. 1 0)-(4. 22).  PLASMA  makes  use  of  two 
further  subprogram  units  PHYCON  and  MAGCON  to  treat  the  case  of  the  plasma 
being  confined  by  the  breech  and  the  case  of  the  plasma  being  confined 
completely  by  the  magnetic  pressure  respectively. 


In  its  present  form,  computer  memory  requirement  by  the  program  is 
minimal  and  is  not  a  consideration  in  implementing  it  on  most  commonly 
available  main  frame  computers.  It  may  also  be  implemented  on  many 
minicomputers  or  microcomputers.  For  some  microcomputers,  memory 
requirements  may  be  a  consideration.  If  so,  Figure  3  provides  a  convenient 
diagram  for  overlaying  the  various  compiled  program  modules  in  linking  them 
together  to  form  a  segmented  running  module. 


6.  CRITIQUE  AND  RECOMMENDATION  FOR  FUTURE  WORK 


As  was  originally  intended,  PARA  was  written  and  used  primarily  for 
planning  and  designing  experiments  aimed  at  elucidating  behavior  of  plasma 
armatures.  It  basically  forms  a  first  generation  code  for  simulating  the 
structure  of  the  plasma  armature.  Though  the  results  of  its  use  have 
exceeded  the  author's  expectations  in  the  extent  by  which  they  agree  with 
experimental  measurements,  qualitatively  and  quantitatively  [15,16],  the  code 
as  a  engineering  design  tool  is  still  deficient  in  many  ways,  both  in  its 
model  for  the  driving  circuit  and  for  the  plasma  armature. 


For  the  driving  circuit,  the  action  of  the  main  switch  and  the  crow¬ 
bar  switch  are  modelled  as  a  fixed  resistance  together  with  fixed  electrode 
potential  drops.  This  is  unsatisfactory  as  the  user  of  the  code  is  required 
to  know  these  values  before  running  the  code.  Further,  the  values  for  these 
parameters  are  likely  to  be  dependent  on  the  current  and  its  time  history. 

They  certainly  depend  on  the  particular  type  of  switch  used.  The  transient 
behavior  of  the  inductor  with  respect  to  its  inductance  is  also  not  modelled. 


For  the  plasma  model,  major  deficiencies  include  the  assumption  of 
mass  constancy,  spatial  uniformity  of  current  density  and  temperature,  and 
single-species  plasma.  None  of  these  assumptions  are  strictly  valid  in 
practice. 


In  railguns,  because  of  the  high  temperature  and  pressure,  ablation 
of  materials  surrounding  the  plasma  may  occur  at  such  a  rate  that  ablation 
products  could  have  significant  effects  on  the  plasma  volume,  temperature, 
resistivity  and  hence  its  resistance. 
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The  role  of  boundary-layer  effects  on  the  plasma  armature  in  the 
railgun  context  are  not  entirely  clear  at  present.  These  effects  are  ignored 
in  the  present  version  of  PARA,  so  also  are  effects  of  higher  dimensions,  2-D 
and  3-D. 

Another  major  deficiency  of  the  code  lies  in  its  use  of  the  concept 
of  electrode  potential  drops  which  are  added  to  the  column  voltage  to  provide 
the  total  voltage  drop  across  the  plasma  armature.  This  is  unsatisfactory 
for  two  reasons: 

(i)  Experimentally,  only  the  total  voltage  drop  across  the  leading  edge 
of  the  plasma  armature  can  be  readily  measured, 

(ii)  Assignment  of  fixed  values  for  the  electrode  potential  drops  in 
running  the  code  can  at  best  be  done  at  present  based  upon 
experience  with  arcs  at  such  lower  currents  (below  ~  10  kA)  and  much 
lower  pressures  (below  ~  1  atmosphere).  Experience  with  these 
relatively  low  current  arcs  show  that  the  30-called  electrode 
potential  drop  varies  between  1  V  and  20  V  depending  upon  the 
conditions  of  the  arcs  [17].  Despite  the  long  existence  of  the 
concept  of  an  electrode  potential  drop,  specific  knowledge  regarding 
its  magnitude,  its  variation  with  current  and  applied  voltage,  and 
its  origin  remain  scarce.  Use  of  the  concept  in  the  present 
version  of  PARA  makes  the  running  of  the  code  somewhat  arbitrary  and 
requires  some  intuition  about  arcs  at  this  level  of  current  density. 


Finally,  by  the  kinetic  assumptions  of  uniform  expansion  and 
contraction  as  well  as  (3.1a)  and  (3.1b),  effects  due  to  MHD  instabilities  are 
ignored. 

In  view  of  these  deficiencies,  the  following  is  recommended: 


(i)  'nie  assumption  of  constant  plasma  mass  should  be  replaced  by  a  model 
for  mass  increase  to  include  ablation  of  rail,  sabot  and  barrel 
materials  due  to  radiative  heating  by  the  arc. 

(ii)  As  railgun  plasma  armatures  are  radiation-dominated  arcs,  more 

accurate  models  for  radiative  heat  transfer  leading  to  non-uniform 
temperature,  pressure,  density,  ionization,  and  current  distribution 
within  the  arc  should  be  developed.  This  may  be  done  in  one  of  two 

ways:  (a)  Approximate  the  process  of  radiative  heat  transfer  within 
the  plasma  as  a  diffusion  process,  with  the  radiative  flux  F_  being 
proportional  to  the  gradient  of  radiation  intensity  I,  that  is 

F  -  XVI 

This  is  the  approach  used  by  Powell  and  Batteh  [9]  and  Powell  [10] 
in  their  1-0  and  2-0  static  analysis  of  the  arc.  (b)  Use  the  exact 
equation  for  radiative  heat  transfer  through  the  use  of  the 
radiation  integral.  The  choice  between  the  two  methods  will  depend 
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on  the  relative  computational  complexity  of  the  two  methods  and  the 
accuracy  of  the  first  method  in  the  railgun  context. 

(iii)  Models  for  the  plasma  sheath  between  the  plasma  column  and  the  rails 
should  be  developed.  Such  models  should  provide  information  on  the 
so-called  electrode  potential  drop,  and  define  more  accurately  the 
transition  between  the  plasma  column  and  the  rails.  In  particular, 
the  transfer  of  heat  between  the  plasma  column  and  the  rails  will  be 
more  accurately  modelled. 

(iv)  The  use  of  more  accurate  expressions  for  plasma  resistivity  taking 
into  consideration  multiple  ionization  and  intense  cross  magnetic 
field  is  warranted. 

(v)  Since  substantial  energy  loss  could  be  caused  by  the  use  of 

inefficient  crow-bar  switches  in  the  case  of  capacitor  bank  driven 
railguns,  provision  should  be  made  in  the  code  for  a  better  model 
for  the  operation  of  the  crow-bar  switch  and  the  coupling  between 
the  crow-bar  circuit  and  the  main  driving  circuit.  Models  for  the 
crow-bar  switch  of  course  depend  on  the  particular  method  of 
switching  and  provision  should  be  made  in  the  form  of  a  user-written 
subroutine.  A  better  model  for  the  main  switch  should  also  be 
used. 

(vi )  For  very  accurate  simulation,  the  transient  behavior  of  the  driving 
Inductor  with  the  consequent  time -dependent  inductance  and 
resistance  should  be  modelled.  This  again  depends  on  the  detailed 
geometry  of  the  inductor. 

(vii)  The  rail  velocity-skin  effect  should  be  modelled.  A  minimum  effort 
here  would  involve  an  integral  model,  at  least  2-D,  which  treats  the 
moving  plasma  column,  plasma  sheath  and  conduction  in  the  rails  as  a 
single  entity  in  a  consistent  way  according  to  the  equations  of 
Maxwell  and  magnetoplasmadynamics. 
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